Vacuum Polarization and Dynamical Chiral Symmetry Breaking: Phase Diagram of 

QED with Four-Fermion Contact Interaction 
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We study chiral symmetry breaking for fundamental charged fermions coupled electromagneti- 
cally to photons with the inclusion of four-fermion contact self-interaction term, characterized by 
coupling strengths a and A, respectively. We employ multiplicatively renormalizable models for the 
photon dressing function and the electron-photon vertex which minimally ensures mass anomalous 
dimension 7, n = 1. Vacuum polarization screens the interaction strength. Consequently, the pattern 
of dynamical mass generation for fermions is characterized by a critical number of massless fermion 
flavors Nf — Nf above which chiral symmetry is restored. This effect is in diametrical opposition to 
the existence of criticality for the minimum interaction strengths, a c and A c , necessary to break chi- 
ral symmetry dynamically. The presence of virtual fermions dictates the nature of phase transition. 
Miransky scaling laws for the electromagnetic interaction strength a and the four-fermion coupling 
A, observed for quenched QED, are replaced by a mean-field power law behavior corresponding to a 
second order phase transition. These results are derived analytically by employing the bifurcation 
analysis, and are later confirmed numerically by solving the original non-linearized gap equation. 
A three dimensional critical surface is drawn in the phase space of (a, A, Nf) to clearly depict the 
interplay of their relative strengths to separate the two phases. We also compute the /3-functions 
((3 a and /?a), and observe that a c and A c are their respective ultraviolet fixed points. The power 
law part of the momentum dependence, describing the mass function, implies f m — 1 + s, which 
reproduces the quenched limit trivially. We also comment on the continuum limit and the triviality 
of QED. 

PACS numbers: 12.20.-m, 11.30.Rd, 11.15.Tk 



Since the works of Maskawa and Nakajima as well as 
the Kiev group Q , it is well known that quenched quan- 
tum electrodynamics (QED) exhibits vacuum rearrange- 
ment, which triggers chiral symmetry breaking when the 
interaction strength a = e 2 /(Air) exceeds a critical value 
a c ~ 1. a c was argued to be an ultraviolet stable fixed 
point defining the continuum limit in supercritical QED. 
Although these works were carried out for the bare ver- 
tex in the Landau gauge, principle qualitative conclusions 
were later shown to be robust even for the most gen- 
eral and sophisticated ansatze put forward henceforth for 
an arbitrary value of the covariant gauge parameter, see 
e.g., Bardeen, Leung and Love |(| demonstrated 

that the composite operator %pip acquires large anoma- 
lous dimensions at a = a c . In fact, the mass anomalous 
dimension was shown to be j m = 1, leading to the fact 
that the four-fermion interaction operator (y^ip) 2 acquires 
the scaling dimension of d = 2(3 — 7 m ) = 4 instead of 6, 
and becomes renormalizable. This is an example of when 
an interaction which is irrelevant in a certain region of 
phase space (perturbative) might become relevant in an- 
other (non perturbative). Consequently, the four-fermion 



contact interaction becomes marginal whose absence can- 
not render QED a closed theory in the strong coupling 
domain. Depending upon the non perturbative details 
of the fermion-boson interaction, it is plausible to have 
7 m > 1, implying d < 4, which would modify the status 
of the four-point operators from marginal to relevant, see, 
e.g., the review article Q, and references therein. The 
upshot of the argument is that the robustness of any 
conclusion about strong QED can be guaranteed only 
if it is supplemented by these perturbatively irrelevant 
operators. Quenched QED with the inclusion of these 
additional operators has been studied in 

Unquenching QED involves inclusion of fermion loops. 
It provides screening and transforms the vacuum char- 
acteristics drastically, changing the Miransky scal- 
ing law for the dynamically generated mass m ~ 
A Exp[— n/^a/a c — 1] to a mean field square-root be- 
havior, i.e., m ~ A ^/a — a c , [HI- See also and 
references therein. Employing a multiplicatively renor- 
malizable photon propagator proposed by Kizilersu and 
Pennington [l^, it has recently been shown that large 
value of Nf restores chiral symmetry above a critical 
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value Nj and the corresponding scaling law itself is a 



square-root: m ~ A ^ 7V| — ATf , [13j. However, these 

results were demonstrated without incorporating four- 
fermion interactions. In this article, we include this ad- 
ditional interaction and establish the robustness of this 
result with the inclusion of all the driving elements which 
influence chiral symmetry breaking, namely, a: the QED 
interaction strength, A: the coupling constant related to 
the four-fermion interactions and Nf. the fermion fla- 
vors whose effect is diametrically opposed to that of a 
and A. We study the details of chiral symmetry break- 
ing in the vicinity of the phase change, mapping out the 
phase space of all these relevant parameters and report 
the results which survive as well as the ones that modify 
in different regimes of this phase transition. 

In section U we introduce the framework of the 
Schwinger-Dyson equations (SDEs), the notation as well 
as the assumptions we employ for our analysis. Section ITT1 
is dedicated to the analytic treatment of the gap equation 
in the neighborhood of the critical plane which separates 
chirally symmetric and asymmetric solutions. Next, in 
section Hill we present the results of our numerical anal- 
ysis. The last section IIVI summarizes our findings and 
provides an outlook for future work. 



I. SDE FOR THE FERMION PROPAGATOR 

The starting point for our analysis is the SDE for the 
electron propagator 



S-\p) = S^\p)+ie 2 



d 4 k 



^s(k)r»(k, P )A^( q ) 



d A k 



Tr[S(k)} 



(1) 



where q = k~p, e is the electromagnetic coupling and Go 
is the four-fermion coupling. We define the dimensionless 
four-fermion coupling A as A/A 2 = Go/(4n 2 ). Sq 1 (p) = 
p 1 is the inverse bare propagator for massless electrons. 
We parameterize the full propagator S(p) in terms of 
the electron wave function renormalization F{p 2 ) and 
the mass function M(p 2 ) as S(p) = F {p 2 ) / (p 1 - M (jp 2 )) . 
^nv(<l) is the full photon propagator which can be con- 
veniently written as 



A M „(g) 



G(q 2 



9^ 



(2) 



where £ is the covariant gauge parameter such that £ = 
corresponds to the Landau gauge. G(q 2 ) is the photon 
renormalization function or the dressing function. The 
full electron photon vertex is represented by T ti (k,p). 
The form of the full vertex is tightly constrained by var- 
ious key properties of the gauge theory, 0, e.g., multi- 
plicative renormalizability of the fermion and the gauge 
boson propagators, [1, Qj|, [HI , perturbation theory, [l5j . 



the requirements of gauge invariance/covariance. 
[20| and, of course, observed phenomenology, [2l|. The 
most general decomposition of this vertex in terms of its 
longitudinal and transverse components is 

4 8 

r*(M = £ + £ ^< (M W> p) , (3) 



where L% = 7", L% = (k + pY{}t + jf), L% = (k + p)^ 
and Li = a» u (k+p) v , where o» v = [7^ ,j v ]/2. The co- 
efficients Xi are determined through the Ward-Takahashi 
identity 



(k- P )^(k,p) = s-Hk) - s^ip) , 



(4) 



relating the electron propagator with the electron-photon 
vertex, [22j. Starting from the limiting form of this iden- 
tity, namely the Ward identity, 



OS-Hp) 

dpn 



r"(p,p) , 



employing the most general form of the fermion propa- 
gator and then generalizing to arbitrarily different mo- 
menta, one obtains 
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F(j> 2 ). 
M(p 2 



F(j> 2 ) 



(6) 



and \^[k,p) = 0. 

It has now been established that the choice of the 
transverse vertex has observable consequences at the 
hadronic level, despite the fact that the simple rainbow- 
ladder truncation is sufficient to reproduce a large body 
of existing experimental data on pseudoscalar and vector 
mesons such as their masses, charge radii, decay con- 
stants and scattering lengths, as well as their form fac- 
tors and the valence quark distribution functions, (23|-|33|. 
For example, the conundrum of mass difference between 
opposite parity states can only be explained through cor- 
rections to the rainbow ladder truncations, (34|, also see 
the review, [35| . In addition to the efforts steered through 
the continuum studies, attempts have also been initiated 
in lattice field theory to compute the transverse form 
factors of the fermion-boson vertex in some simple kine- 
matical regimes, [36l. [37j. Extending these efforts to the 
entire kinematical space of momenta k 2 ,p 2 and q 2 is nu- 
merically challenging and it may require some time be- 
fore the results are made available. However, despite the 
fact that the transverse vertex can have material effect 
on hadronic properties and is crucial in maintaining key 
properties of a quantum field theory, the qualitative be- 
havior of the fermion mass function itself is not signifi- 
cantly sensitive to its details. Therefore, for our purpose, 
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we shall restrict ourselves to the simplest construction 
(Eq. (8) of [l3j|) which, in the quenched limit, renders 
the ultraviolet behavior of M(p 2 ) to be 

M( P 2 ) ~ (^W 2 - 1 ( (7) 

with anomalous mass dimensions j m = 1. This large 
value makes it mandatory to introduce four-point inter- 
actions to ensure self consistency. With this choice of the 
full vertex, we obtain, in the massless limit, 

*v> =(£)"> G (« 2 ) = ($) s > (8) 

where v = a£/(47r) and s = aNf/(3iv). Near criticality, 



where g(p 2 ) = sqG(p 2 ), so = 3a/(47r) and A is the ultra- 
violet cutoff which regularizes the integrals. Note that 
we have employed the simplifying assumption G{q 2 ) = 
G(k 2 ) for k 2 > p 2 and G{q 2 ) = G{p 2 ) for p 2 > k 2 , which 
would allow for the analytic treatment of the linearized 
equation for the mass function as detailed in the following 
section. 



II. ANALYTIC TREATMENT 



Before we venture into the computation of the mass 
function by numerically solving the above non-linear in- 
tegral equation, we find it insightful to make analyti- 
cal inroads. The differential version of the gap equation 
(J9j) simplifies in the neighborhood of the critical coupling 
a c ; viz., the coupling whereat M(p 2 ) ^ solution bi- 
furcates away from the M(p 2 ) = solution, which alone 
is possible in perturbation theory. The behavior of the 
solution near the bifurcation point may be investigated 
by performing functional differentiation of the gap equa- 
tion with respect to M{p 2 ) and evaluating the result at 
Al(p 2 ) = 0. Practically, this amounts to analyzing lin- 
earized form of the original gap equation [i.e., the equa- 
tion obtained by eliminating all terms of quadratic or 



where the generated masses are small, it is reasonable to 
assume that the power law solutions for the propagators 
capture, at least qualitatively, correct description of chi- 
ral symmetry breaking. We choose to study the resulting 
equation for the mass function in the convenient Lan- 
dau gauge. Results for any other gauge can be derived 
by applying the Landau-Khalatnikov-Fradkin transfor- 
mations [20l . l38l . [391 1 or using a vertex ansatz which ef- 
fectively incorporates gauge covariance properties in its 
construction, e.g., 0, 



After taking the trace of Eq. ((T|) , carrying out angular 
integral and Wick rotating to Euclidean space, we obtain 



(9) 

I 

higher order in M(p 2 )]. 

p 2 A 2 

M{p 2 ) = i^ljdk 2 M(k 2 ) + Jdk 2 ^lg(k 2 ) 

p2 

A 2 

+ j^Jdk 2 M{k 2 ) . (10) 


Note that the non-linearized version of this equation 
receives negligible contribution from the region k 2 — > 0, 
while this is not true for Eq. (|10p . This shortcoming 
is readily remedied by introducing an infrared cutoff m 2 
such that M(m) = in. The resultant linearized gap equa- 
tion, Eq. (fT0|) . can now be studied analytically in the 
neighborhood of the critical plane on converting it into a 
second order linear differential equation 

x 2 M"(x) + sxM'(x) + s (l - s)^^ = , (11) 

X s 

with two boundary conditions. Here we have used the 
convenient substitution x = A 2 /p 2 . Following are the 
infrared and ultraviolet boundary conditions 

M'(A 2 /m 2 ) - 0, (12a) 

(i + A)^a. (12b) 

V s / 1-s 

Note that the four-fermion coupling only affects the ul- 
traviolet boundary condition. The differential equation 
itself and the infrared boundary condition do not have 
any direct dependence on it. It is what we expect in- 
tuitively, four-fermion Nambu-Jona-Lasinio (NJL) type 



M{p 2 ) = g( P 2 ) / dk 



M{k 2 ) 
+ M 2 {k 2 



dk 2 



M(k 2 ) 
+ M 2 (k 2 ) 



g(k 2 ) 



A 

A 2 



dk" 



k 2 M(k 2 ) 
k 2 + M 2 {k 2 )' 
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term only generates a constant mass term. Therefore, it 
effectively serves as a cut-off dependent bare mass and 
can hence only influence the dynamics through the ul- 
traviolet boundary condition. 

If we now apply the Lommel transformations: z = Bx 1 
and W(x) = x~ a M(x), the linearized equation can be 
converted into the following Bessel differential equation 

z 2 W"{z) + zW'(z) + {z 2 - A 2 )W{z) = , (13) 

where, 7 = —s/2, a = (1 — s)/2, A = (1 — s)/s, and 
B = y/3a(l — s)/(tts 2 ), where s < 1 . The boundary 
conditions in terms of the function W(z) are given by 

aW{z)+ 1 zW'(z)\ z=B{m/Ay = 0, (14a) 



where Ja(z) and Ya(z) are the Bessel functions of the 
first and the second kind, respectively. The power law 
part of momentum dependence of the mass function 
M{x) = x a W{x) is neatly separated out into the fac- 
tor x a , implying j m = 1 + s. Note that s = 1 corre- 
sponds to a momentum independent photon propagator 
which implies 7 m = 2. Consequently, Eq. ([7|) implies 
that it corresponds to a momentum independent mass 
function, a result which is readily and analytically con- 
firmed from the resultant simple integral equation. This 
is a well known behavior of a contact interaction model 
of the NJL type. Moreover, the quenched limit of 7 m = 1 
is also reproduced trivially. 



a(s ~ X)W(z)- 1 z(s Q + X)W'(z)\ z=B = 0. (14b) 

The general solution of the second order differential equa- 
tion CDS) is 



W(z) = Cl J A (z) + c 2 Y A (z) , 



(15) 



For the homogenous boundary conditions of Eqs. (fT4"|) . 
the non-trivial chirally asymmetric solution of the gap 
equation exists if the following condition holds : 



2aJ A (z) + 72 (Ja-i(z) - Ja+i{z)) 
2aY A (z)+ 1 z (Y A ^(z) - Y A+ i(z)) 



z=B(m/A) s 



(1 + X/s )jB (J A -x{B) - Ja+i{B)) - (1 - s)(l - X/s )Ja(B) 



(1 + X/s ) lB (Ya-i(B) - Y A+1 (B)) - (1 



In the limit of A — > 00, we obtain the following result for the dynamically generated mass m 



S )(l~\/s )Y A (B)- 
(16) 



m 



— = f(a,N f ,X) 



T(^)r(A + 2)2a 

7T7 



(1 + \/s ) 1 B (Ja-i(B) - J A +i(B)) - (1 - s)(l - X/s )J A (B) 



(1 + A/so) lB (Ia-i(-B) - Y A+ i{B)) - (1 - s)(l - X/s )Y A (B) 



.(17) 



Carrying out the Taylor expansion near the critical point, 
we find the following scaling laws : 



m 

A 

m 

A" 

m 

A 



nil . , 

— = A 1 {a-a c 



A 2 {\~\ c ) 



1/2 



1/2 



Ao N 



f 



,V2 



Nf, X fixed (18) 
N f ,a fixed (19) 
A, a fixed (20) 



The momentum dependence of the mass function (based 
upon the numerical calculation discussed in the next sec- 
tion) and the scaling laws for A, a and Nf have been 
plotted in Figs. (flUB]) . 

The critical values of the parameters a, Nf, and A de- 
fine a surface in the 3D phase-space of these parameters. 
The mass function is zero (non-zero) below (above) the 
critical surface, which corresponds to restored (broken) 
chiral symmetry. The analytic expression for the critical 
surface can be obtained by setting mj A = in equation 
(fTT)) . The resultant equation, which can be solved for A, 



is given by 

YiB{J A ^{B) 



A 



■so 



Ja+i(B))-(1-s)J a (B)] 



[ 7 B(Ja-i(S) - Ja+i(B)) + (1 - s)Ja(B)} 



.(21) 



In order to obtain a finite mass for the charged fermion 
in the limit of A — > 00, one requires charge renormal- 
ization. Therefore, in this limit, we impose a(A) = 
a c + to 2 / (A 2 A 2 ). One can thus readily obtain the corre- 
sponding /^-function : 



Pa 



A 



da 
OA 



X,Nt 



-2(a-a c ) 



(22) 



Therefore, a c is the ultraviolet fixed point of fi a , as has 
been observed in [loj]- Identical presence of the fixed 
point for f3\ for A c is readily observed : 



ft 



A 



ox 

dA 



a,N, 



= -2(A - A c ) 



(23) 



The analytical results of this section can be confirmed 
and made precise through a numerical study of the non- 
linearized gap equation ^ . This analysis is presented in 
the next section. 
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III. NUMERICAL RESULTS 



In order to compare and confirm the above analyti- 
cal results, based on the linearized approximation, we 
solve the original non-linear integral equation ([9]) numer- 
ically for varying Nf, a, and A. Depicted in Fig. 1 is 
the fermion mass functions for different values of A for 
a = 2.5 and Nf = 2. The closer we get to the critical 
value A c , the more drastically pronounced is the drop in 
the mass function, indicating the approaching onslaught 
of the phase transition. The scaling law is explored in 
Fig. 2, where the variation of ra/A = M(p 2 = 0)/A with 
A is plotted at the fixed values of a and Nf. The fit of 
the complete numerical data shows that the power of the 
scaling law is slightly different from 0.5. This is expected 
as the mean field scaling behavior (|2"0")) captures the cor- 
rect physics in the immediate vicinity of the critical point , 
only where the linearized version of the equation becomes 
exact. In the same figure, we also superimpose the an- 
alytically derived square-root scaling law which, as ex- 
pected, sits exactly atop the numerical findings in the 
immediate vicinity of the critical point. In Figs. [3] and 
ID we show the variation of the mass section and the cor- 
responding scaling law as a function of a at fixed values 
of Nf = 2 and A = 0.6. These results again confirm 
the validity of the square-root dependence of the dynam- 
ically generated mass on the electromagnetic coupling. 
In Figs. 03 and [6j we plot the mass function in the pres- 
ence of increasing types of virtual fermion— anti-fermion 
pairs and the resulting scaling law as a function of Nf, 
respectively, for a = 2.5 and A = 0.3. These results es- 
tablish the robustness of the conclusions presented in the 
reference [l3T | on the inclusion of the four-point contact 
interaction term. 

Shown in Fig. [7] is the critical curve in A — a plane 
at Nf = 2. The dynamical mass ceases to exist for the 
values of A and a below the curve. The dots in this 
figure represent the points obtained by numerically solv- 
ing the non-linear integral equation of the mass function 
and the solid curve represents the analytical result which 
corresponds to expression (|21[) . This criticality should be 
considered as an extension of its quenched QED counter- 
part obtained in [4l|. For the sake of completeness, in 
Figs. [H] and [51 we present the critical curves in A — Nf 
and a — Nf planes for fixed values of a = 2 and A = 0.3, 
respectively. These figures show that the analytical re- 
sults agree with the numerical findings with a very good 
accuracy. Fig. ^ gives a quantitative picture of how 
the growing number of fermion flavors requires stronger 
electromagnetic coupling. The relation is not linear. The 
screening effect exhausts the strength of the interaction 
faster with increasing Nf. 

Finally, in Fig.[T01 we present the full critical surface in 
the phase space of Nf, a, and A. In the limit of Nf — > 0, 
Miransky scaling law is reproduced. The results pre- 
sented in this paper are qualitatively robust if, instead 
of the multiplicatively renormalizable photon propagator 
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FIG. 1: The mass functions for different values of A at fixed 
a — 2.5 and Nf = 2. Its increasing sensitivity to the variation 
in lambda helps locate the critical strength A c . Dashed and 
solid curves represent mass functions with (vacuum polariza- 
tion of Eq. (|25p) and without feedback (vacuum polarization 
of Eq. ((SJ) from the gap equation, respectively. 



of Eq. (|8]> . we employ any of the following models : 

• One loop perturbative photon propagator, as em- 
ployed in : 



G(q 2 ) 



3tt I A 2 



(24) 



A photon propagator which receives feedback from 
the gap equation, i.e., 



G(q 2 ) = [{q 2 + M 2 (0))/A 



2i* 



(25) 



which is also multiplicatively renormalizable. For 
a comparison, we have also displayed numerical re- 
sults for this latter choice in all the relevant fig- 
ures. As we had anticipated, near criticality, results 
are practically indistinguishable from the ones ob- 
tained from using the model of Eq. ((SJ • 

Note that away from criticality, a complete self con- 
sistent coupled solution of the photon and the fermion 
propagator will be required. However, finiteness of the 
dynamically generated mass for A — > oo forces non per- 
turbative QED to be consistently defined only for those 
values of a and A which lie on the critical surface. This 
is a simple corollary of the argument laid out in @ . Note 
that as we employ a model for the vacuum polarization, 
the running coupling is not our prediction. Following 
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FIG. 2: The scaling law for four-fermion coupling A. Dots, 
solid and dashed lines represent numerical results, fit to the 
numerical data with a power law and analytically predicted 
square-root scaling law, respectively, at a = 2.5 and Nf — 2. 
The mean field behavior of the chiral phase transition is 
evident. Dashed and solid curves represent mass functions 
with (vacuum polarization of Eq. (|25|) 1 and without feed- 
back (vacuum polarization of Eq. ©) from the gap equation, 
respectively. 



Rakow, [42| , if we define the renormalization at q 2 
rather than on the mass shell, we get 



a R {0) = a F%(fl) G R {0) = a G R {0) 



(26) 



because -F'r(O) = 1 for us. Therefore, our model conforms 
to an(0) — > in accordance with the lattice computation 
of unquenched QED, [42|]. As we have argued before, 
when A — > oo, a — > a c . Thus afl(0) = a c Gn(0) — > 0, 
which is associated with the triviality of QED in [42j . We 
use this same model for the vacuum polarization even in 
the presence of the perturbatively irrelevant four-fermion 
interaction terms. This means that on the phase bound- 
ary, the renormalised coupling is zero even in the pres- 
ence of the four-fermion interactions. This is in accor- 
dance with the argument presented in jTTl |. However, 
note that for the practical solution of the gap equation, 
the photon propagator or the running coupling below 
q 2 = K 2 = M 2 (k 2 ) has no bearing on chiral symmetry 
breaking solution. 

We now recall that the Bethe-Salpeter equation gives us 
an approximate relation between the integral over the 
mass function and the "decay constant of the pion (/)" 
given by Eq. (4.42) of [ll| 



[p 2 + M 2 (p 2 )f 



(27) 



FIG. 3: The mass function M(p 2 ) for varying values of the 
electromagnetic coupling a for fixed values of massless fermion 
flavors Nf = 2 and the four-fermion coupling A = 0.6. The 
objective is to hunt a c and determine the nature of the phase 
transition. Dashed and solid curves represent mass functions 
with (vacuum polarization of Eq. (|25|) ~) and without feed- 
back (vacuum polarization of Eq. j8|) from the gap equation, 
respectively. 



We numerically evaluate / for varying values of s and find 
that / — > oo in the limit when the ultraviolet regulator 
A — > oo , a result which suggests that the continuum limit 
is the one of noninteracting bosons, in agreement with 
earlier findings fiol |43| . 



IV. EPILOGUE 

We have studied chiral symmetry breaking for fun- 
damental fermions interacting electromagnetically with 
photons and self interacting through perturbatively 
irrelevant four-fermion contact interaction which is 
required to render QED a closed theory in its strongly 
coupled regime, where, this additional interaction 
becomes marginal and, perhaps, even relevant. We 
have used multiplicatively renormalizable models for 
the photon propagator (with and without the feedback 
from the gap equation) which, we argue, should capture 
the qualitative physics correctly in the vicinity of the 
critical surface in the phase space of (a, A, Nf), marking 
the onslaught of chiral symmetry breaking (restoration). 
The presence of virtual fermion-anti-fermion pairs 
changes the nature of the phase transitions along the a 
and A-axes. The Miransky scaling law softens down to 
a square root mean field scaling behavior as a function 
of all the three parameters a, A and Nf, (|2TJ|) . Study of 
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the mass anomalous dimensions for QED with a model 
vacuum polarization reveals how, quantitatively, the 
momentum dependence of the photon propagator, i.e., 
(p 2 ) 8-1 , filters into the momentum dependence of the 
fermion mass function, namely, (p 2 )( s_1 )/ 2 , through the 
gap equation. We believe that our analysis can and 
should be extended to the study of QCD through its 
SDEs. Situation is ripe for the application of this line of 
approach and technology to QCD, where we are finally 
having the first glimpses of the flavor de pen dence of the 
gluon propagator in the infrared region, [4(|. This is for 
future. 



FIG. 4: The scaling law for the coupling a, investigated 
through the behavior of the mass function near critical- 
ity, Fig. [3j Dots, solid and dashed lines represent numeri- 
cal results, fit of the numerical data to the power law and 
analytically predicted square-root scaling law, respectively at 
Nf — 2 and A = 0.6. Dashed and solid curves represent mass 
functions with (vacuum polarization of Eq. (|25[) ) and with- 
out feedback (vacuum polarization of Eq. ©) from the gap 
equation, respectively. 
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